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Abstract 

The rate of escape of polymers from a two-dimensionally confining potential well has 
been evaluated using self-avoiding as well as ideal chain representations of varying length, 
up to 80 beads. Long timescale Langevin trajectories were calculated using the path 
integral hyperdynamics method to evaluate the escape rate. A minimum is found in 
the rate for self-avoiding polymers of intermediate length while the escape rate decreases 
monotonically with polymer length for ideal polymers. The increase in the rate for long, 
self-avoiding polymers is ascribed to crowding in the potential well which reduces the 
free energy escape barrier. An effective potential curve obtained using the centroid as 
an independent variable was evaluated by thermodynamic averaging and Kramers rate 
theory then applied to estimate the escape rate. While the qualitative features are well 
reproduced by this approach, it significantly overestimates the rate, especially for the 
longer polymers. The reason for this is illustrated by constructing a two-dimensional 
effective energy surface using the radius of gyration as well as the centroid as controlled 
variables. This shows that the description of a transition state dividing surface using only 
the centroid fails to confine the system to the region corresponding to the free energy 
barrier and this problem becomes more pronounced the longer the polymer is. A proper 
definition of a transition state for polymer escape needs to take into account the shape as 
well as the location of the polymer. 


1 Introduction 

Polymer translocation is a common process in various biological systems [T]. A better un- 
derstanding of these processes is important for novel medical applications and treatments as 
well as for new DNA sequencing technology where the molecule is driven through an artifi¬ 
cial channel in a membrane and each nucleotide induces a characteristic currents across the 
membrane from which it can be identified lasiiiiis]. Experimental measurements have also 
shown that it is possible to separate polymers of different lengths and produce an accurate 
drug delivery system based on translocation ilE]. The crossing rate of the polymer can, how¬ 
ever, depend strongly on the environment, thereby affecting the measured signal. A better 
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understanding of the translocation dynamics could help make these types of methods more 
reliable. 

A polymer escaping from a metastable external potential well through a narrow channel 
represents a generic model of such systems. The channel represents both an energetic as well 
as an entropic barrier. The energy barrier can have contributions from steric effects as well as 
the electromagnetic field of the channel and ions in the surrounding liquid. The translocation 
rate of ideal polymers in simple external potentials has been estimated analytically in limits 
where the chain is either significantly larger or smaller than the external potential well Eli. 
An analogy with semiclassical treatment of quantum tunneling of a particle has been used 
m- Other types of polymers, such as ring polymers have also been studied 

For more complex polymers, numerical simulations can provide an estimate of the escape 
rate. But, since polymer escape is typically a rare event on the time scale of atomic vibrations, 
a direct numerical solution to the equations of motion (’molecular dynamics’, MD) becomes 
impractical. The timescale difference can amount to many orders of magnitude. The path 
integral hyperdynamics (PIHD) method [TSl [IS] makes it possible to accelerate the escape by 
applying an artificial bias force, thereby reducing the time interval that needs to be simulated, 
and then subsequently correcting the calculated rate to give an estimate of the true rate. This 
methodology has previously been applied to polymer escape from a one-dimensional external 
potential [TH|. PIHD has been shown to work even for a time dependent bias force [T7]. We 
note that this method is different from Voter’s hyperdynamics method [15], where the bias 
potential is designed to vanish at first order saddle points. 

Alternatively, an estimate based on statistical mechanics rather than dynamical trajecto¬ 
ries started at the initial state can be used to estimate transition rates, if the initial state is 
assumed to reach and maintain equilibrium distribution of energy in all degrees of freedom. 
Such a rate theory approach involves much less computational effort than simulation of trajec¬ 
tories. Kramers theory of chemical reactions in solutions m has, for example, been applied 
to estimate the escape rate of polymers. It assumes, however, a one-dimensional reaction 
coordinate and the question is how to define such a coordinate in a system with multiple, 
coupled degrees of freedom as in a polymer. 

In this article, we present studies of the escape of polymers from a two-dimensional external 
potential well using standard Langevin dynamics, PIHD, and Kramers rate theory. Two 
different polymer models are studied: (1) Self-avoiding polymers without bending stiffness, 
and (2) ideal polymers without excluded volume. The polymers are modeled using the finite 
extension nonlinear elastic (FENE) model with Lennard-Jones repulsive interaction. 

The article is organized as follows: In the following section, the methodology is presented, 
including a description of the polymer models, the PIHD method, Langevin dynamics, and 
the rate theory. In section 3, the results are presented, followed by a discussion in Section 4. 

2 Methodology 

2.1 Polymer models 

The polymers were modeled as strings of beads coupled with an interaction potential and 
subject to a two-dimensional external potential. The Hamiltonian is 


N 



( 1 ) 


i=l 
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where m is the mass of a bead and N is the number of beads in the polymer. The position 
of bead i is given by and the velocity by v^. The interaction potential is 

N 

Vxt(ri) + C/int({rO), (2) 

1=1 


where the external potential is 


^xt (^5 2/) 


+ y‘^), x<xo\ 

Ay — \ojI{x — xt,y + ^coqu"^, X > xq. 


(3) 


The quantities cjq and give the curvature of the well and of the barrier, respectively, 
/Sy the height of the barrier located at and the xq the cross-over point between the two 
parabolas. The external potential is illustrated in Fig. [I](a). A potential function with the 
same x-dependence but without confinement in the ^/-direction was used in the work of Shin 
et al [IH] . 

The interaction between the beads is given by 


N-i N 



= E ^FENEdy - y+1 ) + E ^Lj( y - I'iD) 

(4) 


* <j,i> 


where 

C^fene(^) = ln(l — r^/i?o)5 

(5) 

and 

U'Liix) = 4e[(CT/r)^^ - (a/r)® + 1/4], 

(6) 


The repulsive interaction between the beads is a Lennard-Jones (LJ) potential that is trun¬ 
cated and shifted so that ULj{r) = 0 if r > The shift by e ensures continuity of the 

function. The interaction potential between adjacent beads, which is illustrated in Fig. [I](b), 
also includes an attractive interaction, the so-called finite extension nonlinear elastic (FENE) 
interaction. Non-adjacent beads only repel each other through the Ulj potential, the second 
sum in Eq. (jH) then including all pairs of beads. We will refer to this full interaction model 
as the ’self-avoiding’ polymer. For comparison, we have also carried out simulations with a 
simpler model where non-adjacent beads do not interact at all. The summation over Ulj in 
Eq. (jlj) then includes only adjacent beads. We will refer to this simpler model as the ’ideal 
polymer’. 

2.2 Dynamics 

The dynamics of the polymer is given by the Langevin equation 

mvi{t) + 7ri(t) + Vi$({ri}) = Sj(t), (7) 

where 7 is the friction coefficient, $({rj}) is the total potential energy given by Eq. Q, 
is the gradient taken w.r.t the coordinates of the ith bead and is the velocity of bead 
i. A Gaussian random force H(t) describes the effect of collisions by solvent molecules and 
is defined in such a way that (S(t)) = 0 and = 2^kBT5^^y5{t — t'). Here 
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Figure 1: (a) Illustration of the external potential V^xt(^5 2/); dS])- (b) The interaetion 

potential ?7int between adjaeent beads, see Eq. du. The LJ potential aets between all the 
beads in the self-avoiding ehain model but only between eonseeutive beads in the ideal ehain 
model The C/fene eomponent of the potential diverges at i?o = 2.0; setting a maximum in 
the separation between eonseeutive beads. 


(...) denotes the ensemble average, /i and n Cartesian coordinate indices, ks the Boltzmann 
constant, T the temperature, 6{t) the Dirac delta function and Kronecker’s delta. 

The escape rate is defined as the derivative of the escape probability IZ — dP^^d^) /dt. For 
a single particle, the escape probability can be written using the path integral formulation as 

da [20] 

Pescit) = [ dr f drQP{rQ)P{rQ,tQ\r,t), (8) 

dxf>X}j Jxo<Xl) 

where is the position of the barrier top, P{ro) is the Boltzmann distribution of the initial 
configurations and 

P{ro,to\r,t) J[Dr]exp{-/3I[r{t)]), (9) 

is the probability that the particle has moved from vq at time to to r at time t. J[Dr] refers 
to integration over all possible paths between vq and r, C being the normalization constant 
and I[r{t)] being the action of the path 

I[r{t)] = j dt[mr{t) + ^r{t) + Vi^{r{t))f (10) 

For a polymer, an escape event is defined to have occurred when the x-coordinate of the 
centroid (C) of the polymer, xq = (1/^) has advanced well beyond the location of 

the barrier maximum, beyond x = + 4. Eq. ([HD can be evaluated numerically as 

^ -^traj 

Pesc{t) = "TT ~ U)- (11) 

^traj ^ 

where 0(...) is the Heaviside function, E is the escape time of the ith trajectory and A^traj 
is the total number of simulated trajectories. 
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In PIHD an artificial bias potential Vi){x) is added to the external potential. The action 
can then be split into two parts I[r{t)] = + Is[r{t)]^ where is the action for the 

system in the presence of the bias potential and 

Imt)] = / dt'WUr) • [Vyb(r) - 2B{t% (12) 

is the PIHD correction factor for each trajectory. For the biased system, Eq. m becomes 

^ -^traj 

Pesc{t) = jjz ~ exp(-^/s[ri(t)]), (13) 

"" i=l 

where ti is the escape time of trajectory i [TH] and A^traj is the total number of simulated 
trajectories [20, page 25]. Eq. ffTHl) gives the transition probability for an unbiased sys¬ 
tem in terms of crossing probability obtained from trajectories of the biased system. The 
normalization factor is 

-^traj 

A/'s = exp(-/3/s[ri(t)]). (14) 

i=l 

The bias potential for each bead was chosen here to be 14(x) = when x < x^ and 

Vi){x) = —h/Sy + ^huj‘l{x — xi)Y when x > xq, where 6 is a parameter to be chosen between 
0 < 5 < 1. Thus the bias potential flattens the external potential along x-axis making the 
escape events more frequent. We tried a few different choices of the bias potential, including 
a constant force on all the beads as well as dragging the chain from one end. The one chosen 
here worked best. 

An equilibrium distribution for the initial state, i^(ro), was generated by letting the sys¬ 
tem thermally relax without bias. Configurations were then drawn from this equilibrium 
distribution and the bias potential turned on to generate escape trajectories. Configurations 
were sampled at time intervals of 2r where r is the relaxation time m- 

2.3 Rate theory 

Kramers theory is frequently used to obtain estimates of transition rates for molecules in 
solution [HIES]. It is based on a Langevin description of the dynamics and different expres¬ 
sions for the rate are obtained depending on the magnitude of the friction coefficient. In the 
high friction limit, the Kramers estimate of the crossing rate of a particle escaping from a 
metastable potential is 

Rk = (15) 

27T-fK 

where AK is the height of the energy barrier, ouq is the curvature of the energy surface at the 
initial state minimum and is the magnitude of the negative curvature at the barrier top. 

In order to apply the Kramers formula in the present case, the multiple degrees of freedom 
of the polymer need to be reduced to a single reaction coordinate. One possibility is to choose 
the x-coordinate of the centroid as the independent variable. An effective potential energy 
curve for this one degree of freedom is then obtained by thermally averaging over all the 
other degrees of freedom. The thermal average of a function /(r) for a fixed value of the 
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x-coordinate of the centroid is 


(/)c = 


Zn{xc) 


N N 

/ Yldr'j{{r'i})S{xG - 

^ i=l j=l 




( 16 ) 


where 

N N 

J]^dr'(5(a;c - — (17) 

i=l j=l 

By applying this averaging to the total potential, an effective energy curve $eff(^C 5 ^) = 
($)c is obtained. The friction coefficient for this reduced dimensionality system is 7 eff ~ 

The effective friction coefficient in the Kramers rate expression, 7 ^, was adjusted here to 
obtain a good estimate of the simulated escape rate of a single bead and turned out to be 
7 k = 0 . 827 . A Kramers approximation for the escape rate of a polymer with N beads is thus 
obtained as 

nK{N) = ^-0AE,b(N)_ 

27r7eff 

From the shape of the effective potential curves, the parameters sind 

AF^eff(^) were estimated by fitting parabolas at the initial state minimum and at the barrier 
maximum. 

The internal degrees of freedom of the polymer contribute to the effective potential curve 
$eff- Alternatively, an effective external potential curve without including the interaction 
between beads can be calculated as (xc, N) = (T4xt)c- We compare below the two energy 
curves and the Kramers rate estimates obtained from each one. 

2.4 Simulation parameters 

The values of the various parameters used in the simulations were m = mo = 1870 amu, 
ksT — 1.2 and a — 1.02 nm, which corresponds roughly to three base pairs of DNA. These 
parameters fix the mass, length and energy scales resulting in a time scale characteristic of 
the LJ potential as ^lj = ^/rno^Je = 30.9 ps, where e = 1 ksT. The external potential was 
defined by parameters uq = 0.0014, ui, = 0.032, AV = 0.3 ksT. The barrier was located at 
position xi) — 16 and the crossover between the two parabolas at xo — 12. The parameter in 
the FENE spring constant was kp = 15 and the maximum FENE separation Ro — 2.0. The 
Langevin equation was integrated in time using a velocity Verlet type algorithm [23] which is 
particularly well suited for PIHD. The effective potential curves $eff ^.nd were sampled 
during the escape simulations. The PIHD bias parameter b was chosen between 0.7... 0.9. 

The chosen friction in the Langevin dynamics, 7 = 0.7 (= 3.8 x 10“^ kg/s), corresponds 
to the high friction range of Kramers’ theory 7 ^ m- With this choice of friction, 
viscosity of the fluid surrounding polymer can be estimated to be 7 1.3 x 10 “^ g(cm s)“^ 

(for water 77 = 1 g (cm s)“^) [23]. 

3 Results 

The escape probability was calculated for polymers with up to A = 80 beads using PIHD 
and Eq. m- A linear least squares fit to the calculated Pesc{t) was then used to estimate 


Zn{xc) = J 
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Figure 2: (a) Comparison of the escape rate of self-avoiding polymers from one- and two- 
dimensionally confining potentials. Red line: A potential well with two-dimensional confine¬ 
ment by Eq. Blue line: A potential well with with the same x dependence but no con¬ 

finement in the y-direction (taken from JT^). The added confinement in the two-dimensional 
well results in enhanced escape rates for the longer polymers, (b) Contribution of the external 
potential to the energy barrier in the potential with confinement in y-direction (red lines) and 
without (blue lines). Circles present the effective external potential Kff (xc) in the well bottom 
xq — {) and triangles at the barrier top xq = xi^. 

the escape rate Tl{N). A comparison of the present simulations using a two-dimensionally 
confined external potential and previously reported simulations m for an external potential 
without confinement in ^-direction (similar as Eq. ^ but with no terms depending on y) 
is shown in Fig. [2 The added confinement in the two-dimensional well results in enhanced 
escape rate for the longer polymers and compared to the one-dimensional case. A pronounced 
minimum in the rate is observed at around N = 30. This occurs because the repulsive 
interactions between beads in the longer polymers raise the free energy of the initial state 
with respect to the transition state. 

The total energy depends strongly on chain length with lower bound $eff > 9(A^ — 1) 
but the values of in Fig. [2] reveals that when two-dimensional confinement is present 
Vq — increases faster. This is due to ’’crowding” in the well; the longer self-avoiding 

chains cannot fit into the initial state any more so they experience much higher external 
potential. In the potential without confinement, the energy barrier Kff (x^) starts to decrease 
after N = 32 corresponding the rate minimum which implicates that the chain is more 
elongated over the barrier when ^-directional confinement is not present. 

When the repulsive interaction between non-adjacent beads is turned off (no excluded 
volume interactions), in the ideal polymers, see Fig. El this minimum disappears. In this 
case, the escape rate continues to drop past N = 30. This shows that the reason for the 
minimum observed in the escape rate of the self-avoiding polymers is due to the repulsive 
interactions between non-adjacent beads. The simulations of the ideal polymers were carried 
out up to N = 48, beyond which the free energy barrier becomes so large that even with 
PIHD the simulations become excessively long. 

PIHD is, nevertheless, found to be efficient for the ideal chains. MD simulations with 
2 X 10^ trajectories and chain length N = 32 give a root mean square error of (Jmd = 1.4 x 10“^ 
for the rate, while PIHD with same parameters gives (Jpihd = 0.9 x 10“^. Knowing that root 
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Figure 3: Escape rate for (a) ideal polymers, and (b) self-avoiding polymers calculated using 
PIHD and estimated using Kramers theory, Eq. (HHD. The escape rate of the ideal polymers 
is monotonically decreasing as a function of N up to the maximum length simulated, while 
the escape rate of self-avoiding polymers exhibits a minimum around N = 30. The purple dot 
presents the rate computed using the energy barrier averaged over the tilted line in Eig. It 
shows that taking the shape of the polymer into account in reaction coordinate improves the 
rate given by Kramers theory. 


— 1/2 

mean square error scales as a ^ ^traj estimate that using MD solely would need 

approximately 2.5 times more trajectories for same accuracy. For chain length N = K) this 
ratio is approximately 3 while for = 24 it is approximately unity illustrating that PIHD is 
more beneficial for the longer chains. Typical data for P{t) are shown in Fig. [4] for straight 
MD and PIHD simulations for an ideal chain with N — ?>2. 

The effective potential curves obtained by fixing the x-coordinate of the centroid of the 
polymer and thermally averaging over the positions of beads according to Eq. (USD are 
shown in Fig. [5l The barrier to the escape of the self-avoiding polymer in the $efr effective 
potential first increases with length and then decreases while the location of the maximum 
monotonically shifts towards the initial state minimum. This shift is also seen in the barrier 
of the average external potential V^ff. The statistical sampling is easier for the self-avoiding 
polymer since the effective energy barrier is lower and the direct dynamics sample the saddle 
point region better. For N = 80, the effective barrier for the self-avoiding polymer was too 
small to obtain good statistics. 

An estimate of the escape rate using Kramers rate theory is obtained by fitting the effective 
potential curves with parabolas to extract estimates of the parameters 
and A'\4ff(A^) in the rate expression. 

The rate estimates obtained using Kramers theory Eq. ffTS]) applied to the effective po¬ 
tential curves, are compared with the PIHD simulated results in Fig. El The rate theory 
gives behavior qualitatively similar to the simulations for both polymer models but severely 
overestimates the escape rate of the longer polymers, especially the self-avoiding ones. The 
reason for this is an underestimate of the energy barrier to escape in the effective potential 
curves. The use of the centroid coordinate as a reaction coordinate does not confine the poly¬ 
mers to the transition state region which then leads to an underestimate of the energy barrier. 
This can be seen by evaluating an effective potential function of two variables, the radius of 
gyration, Rg, as well as the centroid. Such two-dimensional effective potential surfaces are 
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Figure 4: Probability that an ideal polymer of length N — 32 has eseaped, P{t), ealeulated 
using either MD or PIHD simulations with 2 x 10^ trajeetories. A linear least squares fit to 
the MD simulation data gives the rate TZ — (6.5±1.4) x 10“^ ^ and a fit to the PIHD simulation 
the rate TZ = (6.5 ± 0.9) x 10“^. 




Figure 5: Effeetive potential energy eurves for self-avoiding polymer eseape from a two- 
dimensional well, using the x—eoordinate of the eentroid as a reaetion eoordinate and ther¬ 
mally integrating over all other degrees of freedom, see Eq. dm), (a) Effeetive external po¬ 
tential (xc, A/^). (b) Effeetive full potential energy $eff(^C5 ^) “ ^eff (0, N). Seeond degree 
polynomial fits are shown in red and green for well minima and maxima, respeetively. The 
$eff(^C 5 ^) curves shift upward with N so 4>(0, A^) is subtraeted for better illustration. 
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Figure 6: Contour graph of the effective potential energy surface, $efr(^C5 ^^5 8); for the 
self-avoiding polymer where xq is the x-coordinate of the centroid and Rg is the radius of 
gyration of a polymer with 8 beads. In this case the energy harrier lies close to vertical line 
corresponding to a fixed xq (see dashed line). A constraint based on the centroid alone can 
then he used to define a good reaction coordinate. 


shown in Fig. E]for self-avoiding polymer with N — 8 and in Fig. Hfor iV = 40. While the 
energy ridge for the shorter polymer is aligned with the xc = 16 vertical line, showing that 
a constraint based on the centroid coordinate alone can confine the system at the barrier, 
the ridge for the longer polymer is significantly tilted with respect to the vertical axis. This 
means that a constraint based only on a fixed value of the x-coordinate of the centroid cannot 
constrain the system in the high barrier region. When the thermal averaging of the other 
degrees of freedom is carried out for the longer polymer and xq — 8, the polymer either has 
rather compact configurations with a small value of the radius of gyration, or a significantly 
larger value. The intermediate values that correspond to the energy ridge are rarely sampled, 
as shown in the inset of Fig. 0 When, however, the effective energy curve is defined by 
averaging along a line that is tilted in the {Rg,xc) plane, so as to lie along the energy ridge, 
the system cannot escape the high energy region and the vicinity of the first order saddle 
point is sampled, giving a larger average energy barrier to escape. 

If the energy is averaged over the tilted line in Fig. 0it is = 4.7 being higher than 
A4> = 3.7 in Fig. [5] where the average is taken over the straight line. Using Eq. fITS]) we can 
obtain the corrected estimate for the rate over such barrier by 7^' (40) = = 

2.2 X 10- < 7^ = 5.1 X 10 ^ which is plotted as a purple dot in Fig. [3l This rate is closer to 

the rate by direct simulations 7^md(40) = 1.42 x 10“^. 

4 Discussion 

The results presented here show how increased confinement of the external potential, going 
from the one-dimensionally to a two-dimensionally confining potential, affects the escape rate 
of the polymers, as shown in Fig. [2l The added confinement lowers the number of possible 
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Figure 7: Contour graph of the effective potential energy surface, $eff ^^5 40)^ for the 
self-avoiding polymer where xq is the x-coordinate of the centroid and Rg is the radius of 
gyration of a polymer with 4 O beads. In this case the energy harrier is titled with respect 
to a line of a fixed xq (dashed line) and a constraint based on the centroid alone does not 
give a good reaction coordinate. Thermal sampling along the dashed red line is dominated 
by configurations that either have larger or smaller Rg than the value at the energy harrier 
(dashed red line in the inset), resulting in an underestimate of the energy harrier. Sampling 
within the tilted, dashed black line line, however, confines the system within the barrier region 
and has maximum density at the first order saddle point (dotted black line in the inset) and 
gives a larger value of the activation energy. This shows that a reaction coordinate using the 
centroid alone will give an underestimate of the activation energy for polymers of this length 
and that a good reaction coordinate needs to be defined both in terms of the location and shape 
of the polymer. 
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configurations that a self-avoiding chain can take when it is sitting in the minimum, increasing 
the free energy of the initial state with respect to the transition state. The curvature of the 
external potential in the added dimension, the ^/-direction, is the same at the barrier and at 
the initial state minimum, and thus this effect is relatively more important in the initial state 
since the polymer tends to be elongated at the barrier. 

A clear minimum in the escape rate of self-avoiding polymers is obtained for intermediate 
length, about N = 30, and this becomes even more pronounced in the two-dimensional case. 
The escape rate of ideal polymers, where repulsive interaction between the non-adjacent 
beads has been turned off, does not show such a minimum for the range studied here. This is 
consistent with the interpretation that the crowding of the beads in the self-avoiding polymers 
in the initial state well is responsible for lowering the free energy barrier for escape. In the 
case of ideal polymers, such crowding effects are largely absent since only adjacent beads are 
subject to a repulsive interaction. 

Qualitatively correct trends are obtained by applying Kramers rate theory to a one¬ 
dimensional reaction coordinate defined as the x-component of the centroid coordinate. How¬ 
ever, the escape rate is overestimated for the longer polymers. For the ideal polymers, Kramers 
rate theory gives closer agreement with the PIHD simulations. This overestimate of the rate 
for the self-avoiding polymers can also be somewhat reduced by defining the effective poten¬ 
tial curve as a thermal average of the external potential only, '\4ff(xc,A^) = (V^xt)c- The 
results are shown in Fig. [3l This has almost no effect on the rate estimated for the ideal poly¬ 
mers, but significantly reduces that of the longer self-avoiding polymers, essentially through 
cancellation of errors. 

The results presented here illustrate that a centroid coordinate cannot give a good reaction 
coordinate for the longer polymers, as also concluded by Debnath and coworkers M- A 
similar problem in defining a one dimensional reaction coordinate has been discussed in the 
context of the polymer reversal problem [SUES]. A good reaction coordinate needs to include 
information about the shape as well as the location of the polymer at the transition state. 
This result is similar to what has been concluded in quantum mechanical rate theory where 
Feynman path integrals are used to represent quantum delocalization. Here, more beads need 
to be introduced in the path integrals the lower the temperature becomes, so an analogy exists 
between reduced temperature in the tunneling problems and length of polymers in the classical 
polymer escape problems. Calculations of tunneling rates using the centroid as reaction 
coordinate have, indeed, shown an unphysical increase in tunneling rate as temperature is 
lowered, see for example Ref. 123 and a good quantum transition state needs to be defined 
in terms of both location and shape of the Feynman paths [2H1- In the present case, a linear 
combination of the centroid coordinate and radius of gyration could be used as a good reaction 
coordinate, but the proper combination of the two will depend on the length of the polymer. 
A systematic optimization of the location and orientation of a hyper-planar dividing surface 
so as to maximize the transition state free energy could possibly be used for the polymer 
escape problem, analogous to what has been done for diffusion problems 123 EOl, It may also 
turn out that a hyper-planar dividing surface does not provide sufficient flexibility to confine 
the longer polymers to the bottleneck region for the escape. This will be studied in future 
work. 
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